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, Abstract. In this paper, we investigate the quantum transport of a double quantum dot coupled with a 

• nanomechanical resonator at arbitrary strong electron-phonon coupling regimes. We employ the generalized 

' quantum master equation to study full counting statistics of currents. We demonstrate the coherent phonon 

^ , states method can be applied to decouple the electron-phonon interaction non-perturbatively. With the 

help of this non-perturbative treatment of electron-phonon couplings, we find that the phonon-assisted 
resonant tunneling emerges when the excess energy from the left quantum dot to the right one can excite 
integer number of phonons and multi-phonon excitations can enhance the transport in strong electron- 
phonon coupling regime. Moreover, we find that as the electron-phonon coupling increases, it first plays a 
constructive role to assist the transport, and then plays the role of scattering and strongly represses the 
, transport. 

Q ' PACS. 73.63.Kv Quantum dots - 71.38.-k Polarons and electron-phonon interactions - 72.70.-|-m Noise 

, processes and phenomena - 73.23.-b Electronic transport in mesoscopic systems 

■ 1 Introduction way, as detectors with high precisions for the quantiza- 
' tion of the resonator positions [10,18,19,20,21]. Therefore, 
, With the increasing promotion of nanotechnology, peo- it is definitely important to unravel the underlying inter- 

00 ' pie now have the abilities to fabricate fertile atomic and play effect between electronics and vibrations on transport 

] molecular quantum devices, such as single superconduct- properties for future quantum technology development, 
ing electron transistors [1], quantum dots [2] and single 

molecular transistors [3]. As a consequence, great atten- In the other hand, the main transport property, cur- 

tions have recently been attracted on research fields of rent, has been extensively studied in the past decades, 

molecular electronics [4.5,6,7] and nanoelcctromechanical However, current noise is not yet fully exploited as a pow- 

systems [8,9,10], because of their immense potentials to erful tool to extract more information to understand how 

bridge scientific researches and industry applications. One the environment affects the transport [22,23]. Especially, 

' of such promising devices is Double Quantum Dot (DQD) at low temperatures, shot noise [24] is the main source 

. system [11], which, as artificial atoms, is able to confine contributing to current fluctuations because of the dis- 

■ one or several electrons to form effective two level sys- creteness of transferring qubits and charges. Full count- 
tem, so as to manipulate and control coherent transport ing statistics (PCS) [25], initially proposed by Levitov et 
of qubits by superposition of individual electrons' natural al. [26] in mesoscopic physics, is proved as a splendid di- 
states [12,13,14,15]. agnostic tool to investigate complete information of quan- 

It was initially observed by Fujisawa et al. experimen- turn fluctuations both in calculations [27,28,29,30,31] and 

tally [16] and then explained by Brandes et al. theoret- measurements [32,33,34]. The main idea of FCS is, by in- 

ically [17] that single electron tunneling on such quan- troducing auxiliary counting parameters, to evaluate the 

tum dots system inevitably interplays with its environ- cumulant generating function of probability distribution 

ment, such as electron-phonon (e-ph) interaction due to of transmitted charges. Current and shot noise are just the 

vibrations of the system. The transport properties are ex- first two cumulants. The higher order cumulants indeed 

tremely sensitive to the motion of the nanomechanical res- show more details to expand our views to study quantum 

onator such that they can be also utilized, in a reverse fiuctuation effects [35,36,37]. The underlying mathemat- 
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ics of FCS is named large deviation theory [38], whieh has 
been widely used in various fields [39,40]. 

In this paper, we focus on effects of quantized vibra- 
tions of the nanomechanical resonator on transport prop- 
erties of the DQD with large voltage bias between two 
electrodes. The vibration is specified by a single phonon 
mode, which interacts with an external thermal phonon 
environment. We derive a generalized quantum master 
equation (QME) including auxiliary counting parameters 
under the Born-Markov approximation and second-order 
perturbation of system-bath couplings. Though two-level 
system strongly coupled to the oscillator was studied in- 
tensively, people usually applied rotating wave approxima- 
tion (RWA) or only investigated the ground state by vari- 
ational methods in the closed systems [41,42]. While ex- 
tended coherent phonon states method [43,44] is adopted 
to capture the characteristics of excitations mixed by elec- 
trons and phonons with arbitrary strong e-ph coupling 
in the absence of RWA. This non-perturbative approach 
renders us beyond the previous perturbation expansion 
method by assuming a weak e-ph coupling [45]. In fact, 
the strong e-ph interaction has already been found in car- 
bon based nanoscale devices, shown in Refs. [46,47,48]. 
In recent quantum electrodynamics experiments, the cou- 
pling of the two-level artificial atom and the resonator also 
reaches the strong regime [49]. Therefore, our work about 
coherent phonon states approach, which is capable of deal- 
ing with arbitrary e-ph coupling strengths, is important 
and favorite. 

The organization of this paper is as follows: we begin 
by modeling the system with generalized QME combined 
with coherent phonon states approach in the framework 
of FCS in section II. Details about how coherent phonon 
states can optimally decouple e-ph interaction non-perturbatively 
will be given there. In section III, results and discussions 
are presented, of which the first three cumulants of prob- 
ability distribution of transmitted charges are scrutinized 
in the full range of e-ph coupling strengths. Finally, sum- 
mary is given in section IV. 




Resonator 



Fig. 1. (Color online) Sketch map of the DQD (among 
three vertical thick lines, specified by two horizontal lines) 
coupled with the nanomechanical resonator (rectangle). 
Left and right areas with slash lines denote electrodes. The 
thermal environment of phonons is depicted by E (circle) . 



is the environment. Vdl describes the coupling between 
DQD and leads, while Vre denotes the coupling between 
resonator and thermal environment. Their detailed expres- 
sions are to be explained in the following. 

We assume only one electron is allowed to stay in DQD 
at most, considering strong coulomb repulsion in quantum 
dots. The system of the DQD coupled with a quantized 
nanomechanical resonator thus is expressed as: 



Hs = 



iz + Aox + go-z{a + a) + cjfjfl a 



(2) 



2 Theory and Methods 

In this section, we first describe the whole Hamiltonian 
and then give the derivation of the generalized QME ac- 
companied by auxiliary counting parameter and FCS. Fi- 
nally, extended coherent phonon states approach will be 
articulated to show its ability to decouple the e-ph inter- 
action. 



2.1 The Whole Hamiltonian 

Hamiltonian of the whole system, as illustrated in Fig. 1, 
is generally described by 



H — Hs + Hi^cad + He + Vi 



DL 



RE, 



(1) 



where Hs denotes the coupled system of DQD and nanome- 
chanical resonator. i?Lcad depicts the two electrodes. H^ 



Oz = \L){L\-\R){R\ and ct^ = + |i?)(L|, where \L) 

(|i?)) is the state of electron occupying left (right) dot. £ 
means the energy level mismatch between quantum dots 
which can be tuned by gate voltage. A is the correspond- 
ing tunneling element between two quantum dots, g is the 
dots-resonator (e-ph) coupling strength. Lo-h denotes the 
quantized mode of the resonator and (a) depicts the 
corresponding phonon creator (annihilator). 

The DQD is connected with two electrodes, described 
by i?Lcad Z]j={L,_R},fc44fe"^J*: '^i*'^ 4 electronic 
energy with momentum fc in j electrode, and c^|j,(cjfc) 
creating (annihilating) the corresponding electrons. The 
Hamiltonian of electrode coupling parts is given as Vdl = 
Ej={L,ij},fc'7fclj)(0|cifc + H.c, where ^{ is the coupling 
strength between DQD and j electrode, and |0) denotes 
the empty electron state of DQD. The resonator is damped 
by external thermal environment, H-^ = X^fc'^fc^I^fc with 
Wfc the phonon frequency and &^(6fc) the corresponding 
creator (annihilator). The damping Hamiltonian is de- 
scribed by Vre — ^k^kd'^bk + H.c, where denotes 
the coupling of resonator and the thermal environment 
with phonon in momentum k. Without loss of generality, 
we set fi,= l,fcB = l,e = l and wt, = 1 as the energy unit. 
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2.2 Generalized QME and FCS 

To measure fluetuations of electron current through the 
right electrode, we add the counting term to the whole 
Hamiltonian as [50]: 



is the contribution from -ffLcad + Vdl , and the contribution 
from He + Vre reads 



= Hs + i?Lcad + He + V^e + VuL (x) , (3) 

where Vol ix) is transformed from the original T^dl , read- 
ing: 



Vt^Ux) = E {e-i>^r^^\R){0\cn,+r^{^\L){0\cLk 

k 



H.c. 



and Npi = c^j.Cijfc is the electron number operator in 
right electrode, x is the auxiliary counting parameter in 
FCS, appearing in the right coupling part of Vdl(x) to 
count the electron number into the right electrode. 

Following the standard procedure, we treat = Vre+ 
Vdl(x) ^ the perturbation. Note here, the couplings be- 
tween the system Hs and its environments are required 
to be weak, but not the e-ph coupling g inside Hs- Ac- 
tually, g can be arbitrary strong, compared to the energy 
scale LOb- In this way, under Born-Markov approximation, 
the modified QME is derived up to the second order of 



dt 



= -z[Hs,p'^it)] 



XJX/E,Load' 



(4) 



for abbreviation, and 



Here Sj — \0){j\ with j g {L, R}, denotes the annihilation 
operator pumping electron on dot j into the j electrode, 
and Sj = \j){0\ is the creator. a+ = a'^{a~ = a) creates 
(annihilates) phonon of the resonator. e„m = £« ~ Cm, 
with e„ describing the nth eigenvalue of Hs and \n) the 
corresponding eigenstate. Wj{e„rn) = -/l,(enm)/j(erim)/2 
and Vjienm) = -f'j(enm)[l - /j(enm)]/2 are the tunnel- 
ing rates of electrons into and out of the DQD, where 

rj{enm) = 2T^J2k hfcP'5(enm " 4) i^ thc spcctral func- 
tion of j electrode and fj depicts Fermi-Dirac distribu- 
tion correspondingly. n{(jJb) = l/[exp(w;,/T) — 1], is the 
Bose-Einstein distribution, with T the temperature of the 
thermal environment. 7b(w) ~ 27r^j, |Afcp(5(aj — Wfc) is 
the spectral function induced by the thermal environment. 
6j^u = 1 ii j = R, otherwise 0. In the following calculation, 
we apply the conventional wide-band limit: 7b(w) = 7b, 
rj{tnm) = rj and consider a large voltage bias to the 
electrodes {pL^e, A^fifi), such that fL^^nm) = 1 and 
fRi^nm) = regardless of the detail information of e„m. 
Furthermore, we set fi(wb) = — 1 — n(— Wf,) = for sim- 
plicity by keeping zero temperature of thermal environ- 
ment, although this constraint can be released. Consider- 
ing = Ij the QME with counting parameters is 
finally obtained by 



dt 



4 



where [A^,B^]^ 

= TrE,Lead(Px) i^ the reduced system density matrix 
composed of DQD and resonator, by tracing out degrees of 
freedom in electrode leads and thermal environment. (0)g ^^^^ 
shows thermal average of O over the thermal reservoir E 
and two electrode leads. The detail derivation of Eq. (4) 
can be found in Appendix A. Note that p^ is a twisted den- 
sity matrix by the counting parameter x- When x = 0, 
reduces to the conventional density matrix. Therefore, for 
the system considered here, we finally obtain the evolution 
equation of p?., as: 



~i[Hs 
2 



(6) 



SlSIp^ - 2SIpISl 



rpI 



"^SiiP^S 



2 ' 



t s 
-a'ap^ 



2apla^ 



p^a^a] 



d Q 

dA 



When X = 0, Eq- (6) reduces to the same equation as 
Eq. (2) in Ref. [45] and Eq. (9) in Ref. [51]. 

Following Ref. [50] , the moment generating function is 
obtained by Q{x,t) = Tr[p^(f;)]. The k-ih order of charge 
fluctuations in the right electrode is derived as {{ANn)^) = 
§-2Ml where ANr = NR{t) - 7Vfl(0). Then by 



-i[Hs,p'^] + L4p^J+La[p'j, 



X=0 



(5) defining thc cumulant generating function of currents: 



where Le[p^] 

E 

j={L,R}r 



V]{.^nm) {^pl\n){n\S]\m){rn\S.j 



H.c. 



H.c 



+ Wi(enm) (^Sj\n){n\S'^^\m){m\p^^ 
^^j(e™n) (\n){n\S]\m){m\plSj +}i.c^ e"^>^*-« 
- i^o{enm) (^,p^|n)(7i|S']|77i)(m| +H.C.) e^-^'^'- 



z{x) 



hm -ln[g(x,i)], 

t— i-oo t 



(7) 



all cumulants of current fiuctuations for FCS can be de- 
duced straightforwardly as 



d{ixY 



(8) 



Fano factor is specified as the ratio of shot noise S'(0)/2 
l'^' and current / = l'-'^': 



m 

2el 



(9) 
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The crucial stop to detect current fluctuations, as we 
shall see, is to derive the cumulant generating function 
shown in Eq. (7). In weak e-ph interactions [45,51], cur- 
rent cumulants have been studied under the Fock space 
of phonons, by considering the perturbation approxima- 
tion of the e-ph coupling strength. For strong e-ph inter- 
action, a huge truncation number of phonons should be 
considered to converge the results in numerics, resulting 
in tough calculations. Hence it is very challenging to deal 
with this system by using the conventional Fock states. 
While the approach of extended coherent states surmounts 
such drawbacks. This method has been successfully ap- 
plied to Dicke model, spin boson model and quantum en- 
tanglement dynamics [43,44]. 

From the definition, the coherent plionon state is com- 
posed by superpositions of infinite Fock states. Therefore, 
the finite coherent states already include infinite phonons 
in Fock space, which makes the coherent phonon basis 
overcomplcte. This overcomplcte property renders us a 
rapid convergence. We apply this method in the follow- 
ing by optimally choosing an effective displacement of the 
resonator compared to the position in the absence of e-ph 
interaction. Then, a new coherent state basis of phonon 
can be constructed to decouple the e-ph interaction. As 
a result, it is expected that coherent states approach is 
more effective than the Fock states one, and can deal with 
arbitrary strong e-ph coupling strengths. 



2.3 Coherent phonon states approach 

We firstly expand the system density operator under elec- 
trons occupation states, and arrange the elements of inter- 
est as one column: p = [poo, Pll. Prr, Prl, Plr)'^ , where 
Pij = («|p^|j>, with belonging to {|0), |L), \R)}. Then 
the evolution is expressed in a clear way 



(10) 























-iA 


lA 








-rii 


iA 


-iA 





-iA 


iA 










iA 


-iA 








where M = P -I- D[-] -I- La[-]. Here, La[-], is the same as 
the last right term of Eq. (6), with replaced by p. For 
P, it is obtained as 



(11) 



which is traditionally used to describe electron transport 
in the absence of e-ph interaction [22]. It is contributed 
by the first three terms in the right hand side of Eq. (6), 
where in the first term —i[Hs,p^, only the pure electron 
part eoz/'i + Agx in Hs is included. 

The contributions of phonon related parts oobO^a + 
gaz{a^ + a) in Hs in the first right term of Eq. (6) leads 
to D[-]. To get the expression of D[-], people usually treat 
the weak e-ph coupling as the perturbation term. How- 
ever, as we will show, if we jump out of the normally used 



Fock states of phonon. and use a modified phonon basis, 
the e-ph coupling can be decoupled non-perturbatively. 
For the case of one electron occupying the L quantum 



dot, LOha^a 



is specified as lo^o^ a + g{a 



a). Then we can define the modified phonon creator and 
annihilator: 



a, 



with a — — , 



(12) 



which are naturally born to decouple the e-ph interac- 
tion term by changing w^a^a -|- g{a^ + a) to the expres- 
sion ojbA''' A — g'^/ujb- ct = g/^b is interpreted as the ef- 
fective displacement deviation of the resonator induced 
by e-ph interaction under specified electron state, com- 
pared to the position in the absence of e-ph coupling. 
Similarly, when one electron occupies the R quantum dot, 
uJba''a + gazia"^ +a) is finally decoupled to utbB^ B — g^ /ub 
by importing — — a and i? = a — a. It is clear that 
with the help of the modified operators having optimally 
chosen displacement a = g/ujb, the decomposition of e- 
ph interaction has been achieved for arbitrary values of 5, 
which is impossible in the conventional Fock space. 

Recall D[-] is derived from a part of —i[Hs,p^, that 
is, from —i[u}ba^ a + gcTzici^ + o)t P^^ when cr^ subjected to 
the electron vacuum state |0), one can easily find JjPoo = 
— iti;b[ata, poo]i similarly as well as for state |L), ^Pll = 
—iijJb[A'^A^PLL\^ and so on. Thus along this direction, we 
obtain D[p] in the new representation: 



D[p] 



-iwftDiag 



7t. 



[a' a, Poo 

[A'^A.pLL 

[B^B,prr] 
B^BpRL ~ PrlA^A 
\A^ApLR- plrB^bJ 



(13) 



where Diag(:) means all elements are in the diagonal po- 
sitions. 

Accordingly, we have electron-states-specified coherent 
phonon basis: 



\n)L 
\n)R 



(at) 



(At)" 



/n! 
(Bt)" 



■|0)o, 
|0)l, 



(14) 
(15) 
(16) 



with |0)i(«) = e~"'/2eT"''>)o, 



where \n)i denotes the nth exciting phonon state, speci- 
fied by the electron state \i) with i = {0,L,i?}. |0)q ^ ^ 

arc the ground states of the resonator depicted by Wfeflta, 
LObA"^ A, LObB^'B for different electrons occupation states. 
It is straightforward to verify that after the effective dis- 
placement shift ±a, A'^{A) with \n)L and B'^{B) with \n)R 
follow the same creation (annihilation) physics and the 
same calculation rules as the conventional un-shifted op- 
erators at(a) with |ri)o- 
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Then, we can insert corresponding coherent phonon 
states to wholly expand p in the form 

P— (poo, Pll, Prr, Prl, Plr)^ , 

with 

Pjz = {■■■ ,i {n\pj,\m)^,---), (17) 

where the index denoting the possible electron state 
{0,L,R} and \n)j,\m)i being the corresponding phonon 
states. Therefore, we can finally derive Eq. (6) completely 
as a master equation: 

|p = Gp, (18) 

where G is the transfer matrix containing counting pa- 
rameters. The details of the matrix elements of the above 
master equation are described in Appendix B. The solu- 
tion of the above evolution equation reads p{t) = e*^*p(0). 
In the long time limit, the evolution is dominated by the 
eigenvalue of G with the largest real part, Ao(x), such that 
p{t) « e^«('')V(0) [31]. Note when x = 0, Ao(x) reduces to 
zero, which corresponds to the steady state of the dynam- 
ics without counting parameters. Therefore, the eigenvalue 
of G with the largest real part is the cumulant generating 
function at the steady state, Z{x) = ^oix)- Then, we can 
apply Eq. (8) and (9) to investigate the quantum trans- 
port properties, by numerical differentiation. 

In previous studies [45,51], current cumulants have been 
derived under bosonic Fock space. However, in that space 
the e-ph interaction could only be tuned in weak e-ph cou- 
pling strength. The main reason is as the e-ph coupling 
strength g is weak, the e-ph interaction term gaz{a^ + a), 
represented in D[-], can be treated as the perturbation. As 
a result, it is straightforward to solve Eq. (6) for the small 
value of e-ph coupling g under Fock space of phonon with 
convergent results by truncating phonon Hilbert space, 
shown in Ref. [45,51]. However, when g is not small, like 
the usual real situations, the perturbation approximation 
of the e-ph interaction is no longer valid. Tremendous high 
phonon states will be excited by the strong e-ph coupling, 
in which case calculations become tough from the perspec- 
tive of Fock states. 

The advantage of coherent phonon states is that the 
modified creator and annihilator favors the displacement 
of resonator induced by e-ph coupling which in turn de- 
couples the e-ph coupling non-pcrturbatively, as we show 
in Eq. (12) and its following discussions. Moreover, from 
Eq. (15) and (16), it is clear that the coherent phonon 
state actually is the superposition of infinite number of 
Fock states, which makes the coherent phonon basis over- 
complete. This overcomplete property renders us rapid 
convergent calculations under the optimally chosen dis- 
placement a = g/ijJb- Therefore, it is natural to choose 
coherent phonon states as the proper basis to investigate 
the effects of e-ph interaction on quantum transports. 

Further, we would like to point out that the extended 
coherent phonon states approach shares the same physics 
with the polaron (canonical, Lang-Firsov) transformation 



■ r'" — \ — -iJl 


g=o.i 
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Fig. 2. (Color online) Fano factor S'(0)/(2e/) versus DQD 
energy gap e with various g. The other parameters are 
given as Z\ = 0.1, Fl = 0.1, Fr = 0.001 and -yt = 0.01. 
a;5 = 1 is set as the energy unit. 

[52] that they both consider the displaced phonon basis. 
While for the detail mathematical treatment, they are 
different. The polaron transformation method is applied 
to decouple the e-ph interaction, but makes the system- 
electrode tunneling terms more complex, by adding a cloud 
of phonons to the operators of electrons, so-called "dressed" 
states [52] . The extended phonon states approach also de- 
couples the e-p coupling; furthermore, it conserves the 
simplicity of the system-electrode tunneling term. Hence, 
it makes the practical calculations of the results (cur- 
rent, cumulants etc.) efRcient and comprehensible. As a 
result, we use the extended coherent phonon states ap- 
proach to investigate the quantum transport properties in 
the present work. 

3 Results and discussions 

In Fig. 2, we first show Fano factor (see Eq. (9)) as com- 
parisons with previous results given in Ref. [45] to vali- 
date our method. When the DQD-resonator coupling is 
weak {g = 0.1,0.2,0.4), we find the results obtained by 
our method are the same as those in Ref. [45]. Moreover, 
our method with coherent phonon states can also work at 
strong e-ph coupling regime (g = 1.0), where the previous 
method with Fock states fails. 

In the following, with the help of coherent phonon 
states, we are capable to scrutinize the DQD's transport 
properties mediated by e-ph coupling in a full range of 
strengths. We focus on analyzing the first three cumulants 
of probability distribution of electron current, though higher 
order cumulants are also available. 

In Fig. 3(a), by tuning the energy level mismatch e and 
e-ph coupling g, we find there exist multiple resonance is- 
lands of current, where l'^^ suddenly becomes large under 
the integer relation e = fewh, fc G N, as is also discovered in 
Refs. [51,53]. This integer relation indicates that the elec- 
tron tunneling is strongly assisted by fc-phonon excitations 
from the mechanical resonator when the energy gap of two 
quantum dots is equal to the energy of k phonons. That 
is, the tunneling from quantum dot L to R will be sig- 
nificantly enhanced if the excess energy can excite integer 
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Fig. 3. (Color online) Current (I'-^-'); Fano factor 
(l(2)/2:(i)) and RTOC are numerically calcu- 

lated in (a), (b) and (c) under parameters of g and e. 
Dashed black line in (b) is used to signify Poissonion trans- 
port with Fano factor equal to 1. The other parameters 
are given as Z\ = 0.1, Fl = 0.1, Fr = 0.001, jb = 0.01 
and cj;, = 1. 



number of phonons of the resonator's vibration, which arc 
then absorbed by the zero temperature thermal environ- 
ment. 

The positive integer resonance is, however, not always 
satisfied in regimes of large energy mismatch of DQD 
(e 3> w&) and weak e-ph coupling (g <^ ujb) where it is 
insufficient to excite extra phonons, or in regimes of ex- 
tremely strong coupling (g 3> uJt) where e-ph coupling 
plays the role of scattering and strongly represses the elec- 
tron current. For negative e, the current is drastically sup- 
pressed by the fact that the resonator has the ability only 
to emit phonons to the thermal environment of zero tem- 
perature. If we increase the temperature of the environ- 
ment, the resonance peak emerging at negative integer 
values of e will be observed. We note that in a different 
setup of triple quantum dots [54], similar integer relation 
for current is exhibited, but as dips (repressions) rather 
than peaks (enhancements). 

Fano factor shown in Fig. 3(b), is more complicated 
than the current. We separate the plot into two sub-regimes 
by the dashed black line, which notifies Poissonion trans- 
port. In the regime surrounded by this line, Fano fac- 




1 

„ 1.4 

-S^ 1.0 



2 4 

e 



Fig. 4. (Color online) Current X^^\ shot noise l'^' and 
Fano factor Z'^^'/l'^^ as a function of e for different g. 
Asymmetric tunneling rates condition: Z\ = 0.1, = 0.1, 
Fji = 0.001 and 7f, = 0.01 for (a), (b) and (c); Symmetric 
tunneling rates condition: A = 0.1, II = 0.01, Fji = 0.01 
and 7f, = 0.05 for (d), (e) and (f). 



tor drops into sub-Poisson regime, where the coupling 
strength g is moderate [45]. The shrink of Fano factor 
mainly results from the dramatic increasing of current in- 
duced by e-ph excitations. When g becomes strong, Fano 
factor rises to super-Poisson regime due to the repression 
of currents. Furthermore, resonance enhancements of Fano 
factor at integer e are clearly observed in both regimes of 
positive and negative e with strong g. 

Fig. 3(c) shows the rcnormalized third order cumulant 
(RTOC), which should not be ignored to signify high or- 
der transport fluctuations [36,37]. For moderate g, RTOC 
is rather small due to resonances of the current. While 
g reaches strong regime, RTOC is strengthened, mainly 
due to the repression of currents shown in Fig. 3(a). The 
multi-phonon-resonance-induced RTOC enhancement still 
appears, which is similar to what happens for shot noise. 

To understand the complicated behaviors of Fano fac- 
tor, we detail the dependence of current fluctuations on 
DQD energy gap e in Fig. 4. Peaks appear at e = kojt 
for both current and shot noise, resulting from the inte- 
ger phonon resonance. When g is moderate (exemplified 
by g = 0.5), oscillatory decays of both current and shot 
noise in Fig. 4(a) (b) and (d)(e) are visible for e > 4.0 
while Fano factor rises gradually from sub-Poisson regime 
to 1.0 with oscillations, see Fig. 4(c) (f). This can be inter- 
preted that for large e compared to g, the e-ph coupling 
is not able to assist the electron transport and tunneling 
events between quantum dots are seldom, so that succes- 
sive transmissions of electrons are almost uncorrelated. As 
a result, the transport dynamics approaches Poissonion 
regime such that S{0) « 2el. When g is strong (exem- 
plified by g = 1.5), current and shot noise are lifted as 
£ increases (see Fig. 4(a) (b)), indicated by the arrows. It 
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and £ = 0, the current can be described [56] by 






(f) 







1.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

g g 



Fig. 5. (Color online) Current l'^^ shot noise I^^-* and 
Fano factor 

as a function of g for some represen- 
tative £. Asymmetric tunneling rates condition: A = 0.1, 
Fl = 0.1, = 0.001 and 76 = 0.01 for (a), (b) and (c); 
Symmetric tunneling rates condition: A = 0.1, /l = 0.01, 
Fr = 0.01 and 76 = 0.05 for (d), (e) and (f). 



suggests that multi-phonon excitations are definitely fa- 
vored in large e-ph coupling g regime to enhance both 
the current and shot noise. Though not depicted here, we 
still observe that current and corresponding shot noise will 
reach the maximum and then decreases when we increase 
e to extremely large values. 

Although the phonon-assisted tunneling that is excited 
by e-ph coupling is not expected at the negative integer 
value of £ at zero temperature, the resonance peaks are 
still observed in Fig. 4(c) (f) for Fano factor. These peaks 
result from the dips of current at the negative integer e, 
shown in Fig. 4(a) (d). The dips are due to the interference 
between the elastic and inelastic transmission scattered 
by the e-ph coupling, instead of opening new tunneling 
channels by e-ph coupling [55]. 

The dependencies of current, shot noise and Fano fac- 
tor on g are detailed in Fig. 5, with four typical transport 
behaviors exemplified therein: 

(1) The first kind behavior appears at e < (exem- 
plified by e ~ —0.5,-1), where current and shot noise 
are both rather small in the whole regime of g. This is 
understandable since in this regime electrons in DQD are 
severely localized, and there is no feedback of phonons 
from noise environment with zero temperature. As a re- 
sult, electrons cannot effectively transfer through DQD by 
phonon-assisted tunneling. 

(2) e ~ depicts the second kind behavior. When 
g < 1.0, current, shot noise and Fano factor are only 
slightly affected, which show that for the DQD with two 
degenerate energy levels, the transport properties are ro- 
bust for weak e-ph couplings. When g increases further, 
the scattering effect of e-ph coupling becomes significant 
and the current decreases. In the absence of e-ph coupling 



Id) 



A'^Fr 



A'^{2 + Fr/Fl) 



(19) 



Consequently, for the case of asymmetric tunneling rates, 
Ki I/j/2 = 0.0005, while for the symmetric tunnel- 
ing case I*-^^ ~ Ff(/3 = 0.0033, consistent with what are 
shown in Fig. 5(a) and (d). In Refs. [57,58], the exact re- 
sult of Fano factor is obtained without considering e-ph 
coupling {g = 0): 



1 - 8A 



.^eHFr/Fl 



1) 



FI^/Fl + SFr/FlA^ 



[F^ + 4£2 



-AA^{Fl/Fr + 2)Y^ 



(20) 

For e = 0, we then can approximately simplify Fano fac- 
tor as 1 for asymmetric tunneling case Fr <ti F^ and as 
I for symmetric tunneling case A ^ Fr = F^. These co- 
incide with the simulation results shown in Fig. 5(c) and 
(f), where Fano factor is nearly 1.0 for the asymmetric 
case, indicating the Poissonion nature while it is strongly 
suppressed around 0.5 when g < 1.0 for the symmetric 
tunneling. 

(3) The third kind behavior is for positive but non- 
integer £ (exemplified by £ = 0.5, 1.5). As shown in Fig. 5(a) 
and (d), these two curves of currents first increase when 
g increases as a consequence of phonon-assisted tunnel- 
ing. An then, when the e-ph coupling strength g is further 
increased, the currents decrease because of the phonon- 
scattering-induced localization. In Fig. 5(b) (e), shot noise 
changes almost synchronously with currents, however, with 
smaller amplitude variations compared to currents. Hence, 
Fano factor experiences from sub-Poisson to Poisson regime, 
gaining a global minimum around 0.5. 

(4) Positive integer e (exemplified by e = 1,3) de- 
picts the forth kind behavior, that is, the resonant elec- 
tron tunneling regime. The e-ph coupling first plays a con- 
structive role to assist electron tunneling. While, as the 
coupling strength g increases, the resonator's vibration 
scatters electrons and suppresses the current dramatically. 
This is similar to the non-integer e cases and the differ- 
ence is that the maximum of the resonant cases is more 
like a platform. Generally, near the resonances, one can 
write [59]: 



Id) = 



de 



FlFr 



27: (£ - fccjb)2 + {Fl + FrY/A 



[/L-/fl].(21) 



m 



This can be successively approximated as I*-^-* ~ rL+Pr,. 
the weak F]^(^rj limit, since the Lorentz function reduces 
to the delta function. Therefore, for asymmetric tunneling 
case Fr <ti Fl, the resonant current nearly reaches the 
maximum l'^^ = Fr and for symmetric tunneling case 
Fr = Fl, the resonant current nearly reaches the maxi- 
mum 

^'■^^ = rR/2, which are consistent with our calcu- 
lations shown in Fig. 5(a) and (d). Furthermore, different 
from the non-integer e cases, for integer e Fano factor ex- 
hibits two local minimums at its minimum basin in the 
sub-Poisson regime, as shown in Fig. 5(c) and (f). 
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Fig. 6. (Color online) {(7z{a^ + o,)) as a function of g and 
e. Other parameters are A ~ 0.1, Fl ~ 0.1, Fji = 0.001, 
7b = 0.01. 

In Fig. 6, we show the expectation of the e-ph inter- 
action term, (crz(a^ + a)), to exam the role of c-ph inter- 
action in the transport. {crz{ci' + a)) = + a)) ^ 
{\R){R\{a^ +a)) describes the discrepancy of resonator dis- 
placement's contributions on electron occupancies on the 
left and right quantum dots. Under resonance conditions 
(e = fcwh), the discrepancy {az{a^ + a)) is rather small, 
compared to off-resonance cases with the same g, which 
indicates the electron occupation probabilities on the two 
quantum dots are almost equivalent. That is, the phonon 
emission of the resonator into the thermal environment 
assists electron transferring from the left quantum dot 
to the right one, rather than blocks electron transport. 
This is consistent with multi-phonon-assisted resonance 
discussed in Fig. 3 and 4. 

However, when the e-ph coupling g becomes strong, 
the discrepancy {az{a^ + becomes significant for given 
£, which implies the phonon assistance mechanism is de- 
pressed. The e-ph interaction no longer assists electron 
tunneling, but plays the role of scattering. Hence, the elec- 
tron is mainly localized in the left quantum dot (indicated 
by the significantly positive value of {az{a'' + a))) due to 
the e-ph scattering, which results in the repression of the 
electron current, as shown in Figs. 3 and 5. 



4 Summary 

To conclude, we have applied the modified Born-Markov 
quantum master equation to investigate full counting statis- 
tics of electron transport through the double quantum 
dots coupled with a nanomechanical resonator. Particu- 
larly, with the help of coherent phonon states method, we 
are capable to non-perturbatively treat the e-ph coupling 
of arbitrary strong strength, with excellent convergence. 
We have shown the first three order current fluctuations in 
detail to signify the non-Poissonion electron transport fea- 
tures in different dynamics regimes. Conditional phonon- 



assisted resonant tunneling has been observed under the 
positive integer relation of e = fcwb, k £ N. We have further 
found that, in strong e-ph coupling regime, multi-phonon 
excitations are favored such that transport is enhanced 
with increasing energy gap e. Moreover, for positive e, we 
have found that as the strength g increases, the e-ph cou- 
pling first plays a constructive role to assist the transport, 
and then plays the role of scattering and strongly represses 
the transport. Finally, we have studied the expectation 
value of DQD-resonator interaction to show the difference 
of contributions between quantum dots, which comple- 
ments out understanding on multi-phonon-assisted elec- 
tron tunneling and phonon-scattering-induced electron lo- 
calization in strong e-ph regimes. 

Though not fully exploited yet, the extended coherent 
phonon states method may be extended to directly in- 
vestigate dissipative DQD system with multiple-phonon- 
electron coupling [44]. We expect our method can be ap- 
plied to other open quantum systems beyond weak e-ph 
interaction constraints, to uncover possible novel many- 
body quantum effects induced by strong e-ph interactions, 
like the avalanche type of transport [60]. 
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A Derivation of Quantum Master Equation 
under Counting Field 

The whole Hamiltonian including system and environ- 
ments under counting field can be described by = 
Ho + Vx given in Eq. (3). Here Hq = Hs + -^Lcad + He 
shows noninteracting term and V-,^ = Vre + Vdl(x) is 
treated as perturbation one, assuming the interacting be- 
tween system and environment/electrodes is weak (not the 
e-ph coupling inside the system). 

Under interacting picture, the motion equation of den- 
sity matrix from the whole model can be shown as 

^ = -^[V,it),piit)l (22) 

where p^(t) = e'"°* p^{t)e-'"°* andV^{t) = e'"°*V^e-'"°* . 
By integrating Eq. (22), density matrix can be expressed 
as 

p^(t) =p^(0)-z f dt'[V^{t'),p'^{t% (23) 
Jo 

= pAO)-^ f dt'[V^{t'),pM 
Jo 

- fdt' f dT[V^{t!),[V^{T),p{{T)]]. 
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Hence the density matrix motion equation in interacting The second block is for the evolution of pll'- 
picture can be derived by 



dp[{t) 



dt 



(24) 



Transforming it back to the Schrodinger picture, one can 
find the reduced system density matrix can be expressed 
by p^{t) = TrE,Lead[e-'^''Vx(^)e'^''*]- As a result, the 
motion equation of the reduced system density matrix is 
shown by 



dt 



-i[Hs,p'^it)] 



dr{[V^,[V^i-T),p^it)] 



xJx/E,Lcad- 



(25) 



dt 



Pn,m = -icJb{n - to)p"„ + Fl ^L{n\k) 00 {l\m)Lpl°i 



k,l 



-tAj2L{n\k)Rpl';-^+tAY,Pn!lR{k\r>i)L 

k k 



- 2v/(n + l)(m+l)p^^i.„+i] 
+njb[-{n + m)p^^^^, + V(" + l)(m + l)p^^i^„+i 

+ Vri^Pn-l,m-l]- 



The third block is for the evolution of pbb- 



Here Born approximation is applied to decompose the 
whole density matrix as Px{t) = p^{t)®PE®PLead - Markov 
approximation should be also considered to replace P^(t) 
in the right second term of Eq. (24) by {t) and extend 
the upper integral limit from t to oo. Besides, [A^^B^-y^ = 
A^B^ — B^A-x and (O)^ Lg^d shows statistical average of 
O over thermal environment and leads. 



d_ 
Jt' 



„RR — A, , („ ^\r,RR r ^RR 



k 



RL A^\^ ^/„\L\ ^ „LR 



Hk)LPk, 



lb. 



+a^p^^t_i - 2V(« + l)(m + l)p,f_fi,„+i] 



B Derivation of Density Matrix Equation 



The fourth block is for the evolution of prl'- 



The evolution Eq. (18) has five blocks specified by the five 
electron states. Let us denote the reduced system density 
matrix elements combined with counting field as 

<™=^■('^l(J|PxN)l"^)^■ 



Then the first block describes the evolution of elements in 
Poo: 



+rnJ2^{^\k)RR{lHoe'''plf 



k.l 



-|[(n + m)pZ. - 2v/(n+l)(m + l)p°%,„+J 



+njb[-{n + m)p°l,^ + v/(n + l)(m + l)p^°+i.„+i 



00 



„RL 



dr 



RR 
Pn.k 



{n\k)Lpt!^m+^b{m-n)p'^t,] 



lb 
2 



-3aVm+ lp^_^+i - as/^p^; 



RL 

?n — 1 



^2v/(n + l)(m+l)pf^i^„+i] 

+n^b[-{n + m. + W)p^i^ 



-2a(V^mp^fi,„ + V^P^.t) 



+2a(Vm + lp^,^+i + xA^p.^,i_i) 

+ v/(n+l)(m + l)p^fi^„+i + V^p^f i^,„_i] 
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So is the last block for plr'- 

d 



k 

rn^m 



k 

lb 



RR 



2 



[{n + m + Aa^)p^% 
ay^p^ 



2 

Li 



n— l.m 



-2V(n + l)(m + l)p^^i^,„+i] 
+n7b[-(?i + m + 4:a'^)p^^„^ 
+2a(^^^p^«l^„, + 
-2a(Vm + + 
+ v/(n + l)(m+l)p^Ji^,„+i + 

For the inner-product of the coherent phonon states in 
these evolution equations, the expression can be deduced 
by following Eq. (14) to (16): 



L{n\m)o 
L{n\m)n 
R{n\m)L 

with Aim (a;) 



o{n\m)R = {-irOnM, 

(-l)'"Am(2a), 

(-l)"A„n(-2a), 



V;^cxp(-y)x"+™ 



min{n,m} 



-2fc 



{n - k)\{m - k)m' 



In the practical calculations, the truncation number 
of phonon occupation states is set as Ntr. Given a spec- 
ified Ntr, there exist 5 (Ntr -1-1)^ equations in all, where 
the 5 characterizes the space dimension of electron states. 
Thanks to the overcomplcte property of coherent phonon 
states, the truncation number of modified phonon space 
has much shrinked even when g becomes large. Through 
our whole work, we set Ntr < 30. We have checked that 
this number is enough and all the results converge with 
relative errors less than 10"^. 
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